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^he  dissociation  of  Bra  in  Ar  was  studied  in  the  same  shoc)c  tube  using 
three  different  time  dependent  observables  which  were  needed  to  measure  the 
rate  of  the  reaction:  (1)  Brj  molecular  absorption,  (2)  Br  atom  two-body 
emission  and  (3)  density  gradient  change,  detected  by  laser  schlieren  teclmique. 
It  was  found  that  the  first  observable  was  most  useful  in  determination  of  dis- 
sociation rate  constants  between  about  1500  and  1800®K  and  the  second  observable 
between  about  1200  and  1700°K.  Use  of  these  two  observables  yielded  rate 
constants  which  were  in  agreement  with  those  earlier  experimental  data,  which 
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are  likely  to  be  most  reliable.  The  laser  schlieren  technique  yielded  new 
dissociation  rate  constants  between  2100  and  3000°K.  These  data  were  found 
to  be  consistent  with  the  "reliable"  "emission"  and  "absorption"  data,  as  well 
as  with  the  earlier  flash  photolysis  data. 

In  order  to  interpret  the  experimental  results  and  to  study  the  dynamics  of 
the  reaction,  the  dissociation  of  Bra  in  Ar  was  studied  by  3-D  classical 
trajectory  calculations  at  1500,  2500,  3500,  and  6000°K.  In  agreement  with 
earlier  trajectory  studies,  it  was  found  that  Bra  molecules  react  only  if 
their  total  energy  is  within  a few  kT  of  dissociation  limit  and  that  metastablc 
molecules,  with  total  energy  above  the  dissociation  limit,  are  particularly 
reactive.  The  average  energy,.  <AE>  , and  the  average  angular  momentum,  <A£>, 
transferred  in  dissociative  and  non-dissociative  collisions  were  calculated 
as  a function  of  total  energy  of  Br2  molecule,  ys^^enerally,  |<AE>|  (non-dissoc  , ) 
and  |<Aii|  (non-dissoc. ) were  found  to  be  considerably  smaller  than  the  same 
quantities  for  dissociative  collisions.  For  energetic  metastable  molecules, 
dissociation  could  be  accomplished  by  a decrease  of  internal  energy  and 
angular  momentum  of  the  molecule  through  collision.  Thus,  collisionally 
induced  rotational  de-excitation  provides  an  additional  mechanism  for 
dissociation  of  diatomic  molecules.  The  non-equilibrium  effects  in 
dissociation  have  been  evaluated  at  3500  and  6000°K  by  the  method  of  multiple 
collisions.  The  calculated  steady  state  dissociation  rate  constants  between 
1500  to  6000°K  are  in  good  agreement  with  the  available  experimental  data 

Finally,  the  validity  of  the  rate  quotient  law,  i.e.  kfQ^-^ard/I^ reverse”  ^g' 
was  demonstrated  by  trajectory  calculation  technique. 
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The  dissociation  of  Br«  in  Ar  was  studied  in  the  same  shock  tube  using  three  dilTerenl  time  de- 
pendent observables  which  were  needed  to  measure  the  rate  of  the  reaction:  (1)  Brj  molecular 
absorption,  (2)  Br  atom  two-body  emission  and  (3)  density  gradient  change,  detecte<i  by  laser 
schlieren  technique.  It  was  found  that  the  first  observable  was  most  useful  in  determination 
of  dissociation  rate  constants  between  about  1500  and  1800°K  and  the  second  observable 
between  about  1200  and  1700°K.  Use  of  these  two  observables  yielded  rate  constants  which 
were  in  agreement  with  those  earlier  experimental  data,  which  are  likely  to  be  most  reliable. 
The  laser  schlieren  technique  yielded  new  dissociation  rate  constants  between  2100  and  3000°K. 
These  data  were  found  to  be  consistent  with  the  “reliable”  “emission”  and  “absorption”  data, 
as  well  as  with  the  earlier  flash  photolysis  data. 

In  order  to  interpret  the  experimental  results  and  to  study  the  dynamics  of  the  reaction,  the 
dissociation  of  Bri  in  Ar  was  studied  by  3-D  clas-sical  trajectory  calculations  at  1500,  2500, 
3500,  and  6000°K.  In  agreement  with  earlier  trajectory  studies,  it  was  found  that  Brj  molecules 
result  only  if  their  total  energy  is  within  a few  kT  of  dissociation  limit  and  that  metastable 
molecules,  with  total  energy  above  the  dissociation  limit,  are  particularly  reactive.  The  average 
energy,  <a£  >,  and  the  average  angular  momentum,  <Af  >,  transferred  in  dissociative  and  non- 
dissociative  collisions  were  calculated  as  a function  of  total  energy  of  Brj  molecule.  Generally, 
|(AE>|  (non-dissoc.)  and  |(Af)|  (non-dissoc.)  were  found  to  be  considerably  smaller  than  the 
same  quantities  for  dissociative  collisions.  For  energetic  metastable  molecules,  dissociation 
could  be  accomplished  by  a decrease  of  internal  energy  and  angidar  momentum  of  the  molecule 
through  collision.  Thus,  collisionally  induced  rotational  de-excitation  provides  an  additional 
mechanism  for  dissociation  of  diatomic  molecules.  The  non-equilibrium  effects  in  di.ssociation 
have  been  evaluated  at  3500  and  OOOC’K  by  the  method  of  multiple  collisions.  The  calculated 
steady  state  dissociation  rate  constants  between  1500  and  6000°K  are  in  goo<i  agreement  with 
the  available  experimental  data. 

Finally,  the  validity  of  the  rate  quotient  law,  i.e.  ktorw«rd/kreve«e  = A'„,  was  demonstrated 
by  trajectory  calculation  technique. 


1.  Introduction 

Various  earlier  investigators  have  demonstrated 
that  the  absolute  magnitude  of  atomic  recom- 
bination rate  constants,  k„  are  unusually  high, 


and  that  they  have  pronounced  negative  tem- 
Iterature  coefficients.  Furthermore,  for  chlorine, 
bromine  and  iodine,  k,,  deduced  from  shock 
wave  dissociation  studies,  were  often  found  to  be 


• Research  sponsored,  in  part,  by  the  Air  Force  Office  of  Scientific  Research,  Air  Force  iSystems 
Command,  USAF,  under  grant  No.  AFOSR-71-2016.  The  United  States  Government  is  entitled  to 
reproduce  and  distribute  reprints  for  Government  purposes,  notwithstanding  any  copyright  notation 
hereon. 
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Fig.  I.  E.xperimeiital  .shock  wave  profiles,  neiierated  on  a computer  usinft  rate  constants  available  in  the 
literature.  /» is  the  photomultiplier  signal  before  Brj  is  ititrodiiced  in  the  shock  tube.  The  initial  temperature 
was  12G4°K. 


higher  than  A,,  measured  by  other  teelmiiiues, 
such  as  flash  photoly.si.s  which  is  useful  in  re- 
combination rate  measurements.  These  latter 
findings  imply  that  the  rate  quotient  law, 
A(forwardi 'fc(reverse)=  A'(eq),  is  violated  for 
the  speeial  case  of  recombination-rlissociation 
reactions.  Such  experimental  data  arc  particularly 
abundant'”'"  for  the  case  of  brondne  reacting 
in  an  excess  of  an  inert  gas. 

The  purjiose  of  the  jiresent  investigation  was 
to  study  the  dissociation  of  bromine  in  argon  in 
order  to 

(1)  extend  dissociation  rate  measurements 
over  a wider  temperature  range  using  shock 
waves, 

(2)  intei|)ret  the  experimental  fimlings  in 
terms  of  elassieal  trajeetory  calculations  and 

(3)  imciiver  possible  caii.ses  for  the  above 
discrepancies. 

2.  Shock  Wave  Study  of  Bri  Dissociation 
Rate  in  Ar 

The  iliscreiiancies  among  various  shock 
wave'''-'’’  and  flash  photolysis"’  data  for  Hr^ 


-b.\r=21b  + .\r  reaction  arc  particularly  pro- 
nounced in  the  neighbourhood  of  1200°K. 
In  order  to  better  understand  the  nature  of  the 
discrepancies,  we  have  generated  on  the  comimter 
the  .shook  wave  jirofiles  at  12f)4°K,  using  the 
experimental  conditions  given  in  Table  III, 
experiment  #8,  ref.  7.  Details  of  such  com- 
liutations  were  discussed  iircviously.'"  The 
experimental  observable  chosen  in  these  com- 
])utations  was  the  optical  absorption,  I,,  and  its 
time  dependence  during  the  reaction.  (Optical 
absorption  was  the  commonly  used  observable 
in  many  experiments.)®”' 

In  Fig.  1,  the  curve  set  1,  upper  line,  rcprc.sents 
the  ideal  base  line,  i.e.  no  reaction.  The  middle 
curve  is  calculated  using  the  rate  constant 
obtained  from  flash  photolysis  experiments"' 
and  ideal  shock  wave  eiiuations.  The  lower 
curve,  set  1,  is  derived  on  the  same  assumptions, 
except  that  the  rate  constant  reported  by 
Warshay'  is  used.  Crosses  would  have  yieldeii 
the  middle  curve,  except  that  they  were  obtained 
by  including  the  boundary  layer  effects.®  Cal- 
culations described  iiy  ojien  triangles  are  the 
same  as  those  deserilH'd  by  crosses,  except  that 
shock  attenuation  of  one  percent  i)cr  meter  is 
as.sumed.  The  latter  two  sets  of  data  indicate 
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t\iat  till'  (liMTopamy  lictwivii  tiii'  miiiillt'  ami 
till'  liiwiT  l urvi's  caiiiiot  lit'  attriliutnl  tii  iKnimlarv 
lai't'r  t'lTt'rl''  ami  l*»  tlii'  >ln>fk  atti'iiiiatinii. 

(’urve  si't  11  nil  l it:.  1 i"  till'  saiiii'  as  tlir  iiiiclillr 
ami  lowiT  curvi's,  <,.(  1^  I'xri'pt  that  O.ri'','  slint 
linisi'’  ill  (/,  /o)  is  as>iiiiu'(l.  ('iirvt'  >cl  III  is  tlif 
saiiii'  as  i-urvi'  snt  11,  I'Xivpt  that  llin  tiiiif  srali' 
is  twicf  as  fast  as  in  set  II.  In  niiliT  In  pri'.'i'iit 
till'  ii'siilts  nil  till'  balin'  urapli,  riirvt"  11  ami  111 
wore  ailiiti'arily  (li'platoil  frniii  tlio  nri(iiii.  wliiili 
aotilally  is  the  nriKin  nf  oiii  ve  sot  1. 

A ooiiiparisnii  nf  tlio  twn  oiirvos,  sot  II,  I'itt.  1, 
sugROsts  that  it  shniilil  have  Ihx'II  pn.ssihlo  tn 
doteot  a ohaiigo  in  rato  liy  a faotnr  nf  throo, 
which  is  tho  disorojiaiioy  Ik'Iwooii  shuck  wave'” 
and  flash  phntnlysis”'  data,  especially  if  the  Iniig 
time  recnrds  were  tn  lie  examined.  Ilnwover,  it 
call  ho  alsn  soon  frniii  I'ig.  I that  s|npc~  frnni 
which  tho  rato  cmistaiits  had  tn  ho  calculatod  in 
tho  oxporinioiits-”  wore  .small.  Tho  haso  lino  in 
such  experiments  could  ho  ill-dofinod  In'cailse 
nf  |inssil)lo  iinporfoctiniis  in  flow  liohind  tho  shuck 
front.  Tho  existence  of  ralidnm  shot  nniso^  ^ 
.shnuld  iiitriMlnco  further  unoortainty  in  rate 
iiioasuronionts.  I'nr  tlio.-o  roasniis,  the  “ahsnrp- 
tinii"  shock  wave  data  appear  tn  In-  relatively  un- 
roliahlo  arnund  120()“K.  Such  arguments  are  not 
valid  if  slopes  nf  tho  ahsnrptinn  iirnfilos  are 
approciahlo  and  can  ho  nioasurod  imanihignnusly . 
Tho  latter  situation  arises  at  higher  toiiipora- 
tnros  ( l.'iOt)  to  1MH)°K),  and  it  is  in  this  range  that 
tho  rate  cnnstaiits,  iihtainod  in  ahsnrptinn 
spoctrn.sonpy  oxporinioiits^  ' ’ are  in  a fair  agroo- 
mont  with  each  other  and  are  eniisi>tent  with 
flash  photolysis  data.'" 

If  the  nh.-orvahle,  used  for  dotormiliation  nf 
rato  constants  in  shock  wave,  is  tho  two-hmly 
omi.s:-ion,' * then  tho  hase  lino,  which  corro- 
spniids  to  zero  omissinii,  is  practically  iinisoloss 
and  is  unamhigunusly  dofinod  in  cxporiiiionts. 
Fiirthormnro,  in  this  ca.so,  tho  .shock  front 
position  is  easily  dotonniiiahlo.  Thoroforo,  tho 
uncertainty  in  tho  rato  is  largely  dotorniinod 
hy  tho  noise  in  tho  signal  itsi'lf,  prnvidisl,  nf 
course,  that  tho  rato  is  not  ton  slow  nr  ton  fast. 
For  those  reasons,  tho  twn-hndy  omissinii  toch- 
liiijilo' ■ should  provide  accurate  rato  cnnstaiits 
at  Inwo.'t  toniporaturos.  In  order  tn  dotormilio 
tho  advantages  and  limitatinns  of  this  "omi.-sion  ' 
toclmiipio,  wo  have  gonoralod  on  the  cninpiitor, 
protilo^  similar  to  that  shown  on  Fig.  1,  except 
that  wo  wore  using  tho  two-l««ly  oiiiis.'inii  as 
an  olisorvahlo.  TIiom'  profiles  suggest  that  tho 
two-lnsly  omis.-inii  toclmiipio  ic  mdoi'il  roliahlo 
at  Inwo't  toniporaturos,  i.o.  at  almut  12IHI°K, 
provided  that  shot  noise  is  not  too  largo.  However, 
It  proved  to  !«■  difliciilt  to  dofiiio  iiiioipiiviM  ally 
tho  initial  rates  at  highe't  toniporaturos,  la'caii'i' 
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the  omission  is  tho  <|iiadratic  fniiclimi  nf  Hr  almii 
concont  ration.  Therefore,  tho  iiio't  roliahlo 
region  for  rato  nioa'iiroinont'  u-iiig  twn-lindy 
oinis>inn  toclmiipio  appears  to  1k'  hotwooii  ap- 
prnxiniatoly  12(H)  and  17(HI°K. 

.Mthniigh  tho  tomporatiiro  range  for  roliahlo 
Iiioasuronionts,  Imtli  hy  "ah>nrptinn  " and 
“oinissinii''  toclmiipios,  might  have  Ik'oii  ang- 
nioiitod  hy  u>o  of  an  iniprovod  apparatus,  hy  a 
larger  shock  tuls',  and  hy  varying  Itr-ilihiont 
ratio  and  Hr-.;  enneontratinn,  it  proved  diliiciilt 
ill  tho  past  to  ohtain  roliahlo  rates  nf  dissociation 
holow  12(M)°K  and  ahovo  alMiiit  2(KHl"K. 

It  is  important  to  note  that  tho  '‘roliahlo  ' 
.shock  Wave  ‘'omission  " and  "ahsorption  " moas- 
uromonts  are  in  a goml  agrooniont  with  fla'h 
photolysis  data.'"  In  order  to  confirm  those 
conclusions,  wo  have  ropoatod  hoth  tho  ‘■omis- 
sion" and  tho  ‘‘ahsorption"  inoasiiromont' 
Ik'Iow  1S(K)'’K,  Using  tho  shock  tilho,  do'crihod 
previously.'  .\s  aiiticipatod,  the  rosiilt-,  not 
shown  hero,  agrood  with  work  of  provioii'  in- 
vosf  igatoi’s. 

I'ndor  those  circimistaiicos,  it  ~oiniod  ap- 
proprialo  to  oxtond  tho  tompoiaturo  range, 
oniploying  the  same  shock  tiiU''  hut  ii'iiig  a 
now  lasi'i-  sohlieron"'  detect  ion  toclmiipio.  It  was 
felt  that  such  an  approacli  coulil  Iwtfcr  ili'finc 
tho  dissociation-roconihination  rato  constant 
over  a wide  toniporaturc  range,  and  would  make  a 
oomparison  hotwi'on  sli,,ck  wave  and  flash 
Iihotolvsis  data  more  inoaiiingfiil. 

In  laser  schlioron  toclmiipio,  tho  post-'liock 
di'iisit.x'  gradii'iit  causi'd  hy  tho  diss(M*iat ion  is 
folloHoil  hy  nioa'iiring  the  dolloctioii  of  a laser 
Iwaiii  as  a function  of  tinio.  In  older  to  comliicl 
a lasi'i'  schlioron  oxpoiiinont,  tho  'hock  tiilie, 
doscrilH'd  previously,'  wa-  fitted  with  a poli'liod 
iiickol-plateil  stix'l  pipe  olio  iiiotor  long  and  throe 
iiiclios  1.1).,  and  a stahlo  Ilo-.No  ( \\  laser  was 
il'od  to  moasiiro  tho  rate.  I'ho  oxpoi  imontal 
si't-iip  was  similar  to  that  of  tho  previous  work- 
(,j.^i;i-K*  .,tid  ^vill  not  ho  doscrihi'd  hero.  In  order  to 
iihtain  rato  constants  at  highest  toniporaturos,  a 
ililiilo  Ci.7'  i ) sohitiiin  of  lir-i  in  .\r  wa-  ii.'od,  Tho 
late  constants,  dcdiicoil  from  tho  raw  data  in  a 
oonvontional  inaimoi " are  shown  in  !■  ig.  2. 

It  can  lie  'oon  that  tho  comhiiiat ion  of  ‘‘roli- 
ahlo" "oiiiis'ion"  data'  iK'twoon  l'2IHIand  ITIHI^K, 
"roliahlo  "ahsorption"  data  hotwia'ii  hilM) 
and  l.‘'l(H)°K  and  now  hi'or  schlioron  data  Iwlwooii 
2KMI  and  ;i(HH)°K  yield'  a smooth  lino  which  is 
consi'tont  with  earlier  Ihi'li  photol\si'  data.'" 
It  can  also  lie  si'oii  (I'ig.  2)  that  tho  pro'ont 
results  suggest  that,  at  'hock  wa%o  tomporatiiro', 
tho  negative  tomporatiiro  cnofliciont  for  tho 
2llr  + .\i.  ‘Hrs-+-.\r  reaction  is  nearly  unit.\, 
which  Is  smaller  than  that  roportod  hy  iireviotl' 
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FiCi.  2.  A comparisoM  of  calculated  and  experimental  recombination  rate  eonstanbs  for  Hr!+Ar  = 2Hr+Ar 
reaction,  lief,  la  presents  .some  additional  twenty  rate  constants,  obtained  by  two-body  emiRsion  technique 
below  1700°K.  which  are  in  anreement  with  Lli  data.  Open  triangles — trtijectory  <'aleiilations  of  recombina- 
tion rate  constants  CRef.  20). 


investigators,’®-’  who  .sugge.sted  tliat  the  nega- 
tive temperattire  coefficient  is  of  the  order  of 
three.''' 


3.  Dynamics  of  Dissociation  Using 
Classical  Trajectory  Method 

.3-1)  classical  trajectory  study  of  Hr  atom  re- 
eoinliination  in  argon  have  already  provided 
useful  information  on  the  dynamics  and  rate 
eonstaiits  of  the  reaction.'’”®  The,  .satne  princi- 
ples.'’”’" cotipled  with  .stratified  and  im|)ortance 
satn|)ling  techni(|ues'"---’  which  improve  the 
statistics  of  calculations,  are  used  here  to  study 
the  dynamics  of  Hi'j  dissociation  in  .Vr. 

In  general,  the  dissociation  rate  coefficient 
at  any  time  t can  be  written  as 

TT/a)"’  / <f>Ui)S{ii)(lii  (1) 

Jk 

where  S(u)  is  the  reaction  cross-section,  which 


is  a function  of  total  internal,  i.e.  vibrational 
plus  rotational,  energy  u.  Here  it  is  averaged 
over  the  relative  translational  energy  at  tempera- 
ture T,  and  <t>  is  the  normalizeil  di.stribution 
function  of  the  reagent,  i.e.  Hrj,  internal  energy: 


In  the  eq.  (1)  the  integration  is  performed 
ovei-  the  entire  internal  energy  .siiaee  /f,  com- 
jiatible  with  a classicall.v  bound  moleculm 
It  i.s  known  from  earlier  studies,'’  that  the 
reagent  molecules  are  reactive  only  if  « exceeds 
a certain  minimum  value.  Therefore,  the  eq. 
(1)  may  be  rewritten  as 


kj=  {iikT  TtiY'-  j ^{u)S(u)  ihi  (2) 

where  the  subscript  /)  indicates  that  the  inte- 
gration is  executed  over  the  entire  reactive 
energy  range.  In  order  to  ealculate  the  rate 
coefficient  by  sampling  only  from  the  reactive 
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Flo.  3.  I)is.sociation  cross-section  of  Brj  in  Ar,  Sd,  as  a function  of  Brj  internal  energy.  Dashed  lines  with 
triangles,  1500°K;  solid  line  with  circles,  2500°K ; dashed  line  with  squares,  3500°K;  solid  line  with  diamonds, 
6000°K. 


energy  range,  the  eq.  (2)  may  be  replaced  by 
kd=  [SkT  <t>(it)S(u)  du  j <t>(u) 


• / 0(w)  du 
Ji> 


= {v)-{S)-  U(u)du  (3) 

Jn 

where  (v)  is  the  average  relative  velocity  and 


(‘S)  is  the  average  dissociation  cross-section 
within  the  entire  reactive  energy  region. 

In  the  subsequent  three  seetions,  the  dynamical 
l)roperties  of  the  HI  IS,  eq.  (3),  arc  investigated 
liy  trajectory  method.  In  Section  7,  the  dissocia- 
tion rate  constants,  calculated  from  the  eq.  (3), 
are  conqiarcd  with  the  experimental  results. 

The  interaction  potential  for  the  system 
.Ar-Br-Br  is  identical  to  that  used  earlier.™  The 
well  depth  and  inter-nuclear  separation  de- 
scribing the  .Ar-Br  interaction,™  i.e.  and 

were  taken  to  be  0..A  Kcal  hiole  and 
3.32  A,  re.speetively.  Five  strata,  1 kT  each, 
between  —3  kT  and  -|-2  kT  around  the  dissocia- 
tion limit,  were  used. 
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u/RT  u/RT 

Fig.  4.  Contribution  to  as  a function  of  internal  energy  of  the  reagent  molecule.  Note  that  in  I (a)  tlie 
histogram  labelled  CK,  enlarged  by  a factor  of  5,  corresponds  to  the  collisional  release  mechanism. 

4.  Effect  of  Bf2  Internal  Energy  on  Dissociation  The  contribution  to  A/‘>,  the  equilibrium  rate 

coefficient"'*  (one  way  flux),  from  various  internal 
The  (li.ssoeiation  eross-seetions  of  Hro  in  Ar,  Sd,  energy  strata  at  different  temperatures  is  .shown 

are  .shown  in  Fig.  3 at  different  internal  energy  in  Fig.  4.  The  A'/'>  may  be  defined  by  e<|.  (3), 

strata  and  at  four  different  temperatures.  The  using  the  assumption  that  the  integral  J<f>(u)  du 

curves  show  that  oidy  the  molecules  close  to  can  be  calculated  from  the  ifoltzman  distrilni- 

the  dissociation  limit  have  an  appreciable  dis-  tion  function  and  that  (.S')=  ).  At  each 

sociation  cro.ss-section.  This  finding  is  in  agree-  temperature,  the  reagent  molecules  with  internal 

ment  with  previous  .study,”  which  indicates  energy  less  than  3kT  below  the  dissociation  limit 

that  nascent  molecules,  formed  in  recombina-  did  not  contribute  to  the  rate.  Hat  her,  the  ludk 

tion,  pos.se.ss  energies  close  to  di.ssociation.  of  di.s.soeiation  involved  inolei'nles  within  ± IkT 

There  is  a large  increase  of  cro.ss-scction,  as  the  of  the  dissociation  limit.  This  result  juslities  Use 

internal  energy  of  Hr^  increases  from  bound  of  e<p  C2),  rather  than  eipiation  (1),  in  our  cal- 

(«<0)  to  mctastable”  (u>0)  energy  region.  culations.  It  also  can  be  seen  from  Fig.  4 that  the 

The  effect  is  especially  jn'oivounced  at  high  tern-  mctastable  Hi'j  molecules  are  just  as  im|iortant 

peratnre.  This  occurs  because  mctasuible  mole-  as  the  bound  Hr^  molecules  in  the  dissociation 

cules  have  a higher  spatial  extension,  md  becau.se  |)rocesses.  Therefore,  one  may  conclude  from  the 

they  dissociate  via  the  rotational  de-excitation  principle  of  detailed  balancing,  that  the  nascent 

mechanism,  to  be  discus.sed  in  section,  ,'),  products  of  recombination  luedominantly  po]m- 
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u/RT 


Fig.  5.  Average  energy  transfer  vs.  reagent  internal  energy.  Dissociative  collisions  are  described  by  filled 
symbols:  heavy  solid  line  with  diamonds,  ()000°K;  heavy  dashes  with  squares,  3500°K;  fine  solid  line  with 
circles,  2500°K;  fine  dashes  with  triangles,  1500°K.  Non-dissociative  collisions  are  described  by  the  corre- 
sponding unfilled  .symbols,  with  no  lines  connecting  them. 


late  the  highly  energetic  states,  of  which  a sub- 
stantial fraction  are  the  nietastable  states.  This 
conclusion  is  consistent  with  the  results  of  pre- 
vious recombination  studies.'’""*  .At  the  same 
time,  this  conclusion  casts  doubt  on  tliose  tra- 
jectory calculations”  which  are  based  on  the 
a.s.sumi)tion  that  deactivation  of  metastable 
molecules  determines  the  rate  of  recombination. 

In  the  present  calculations,  the  collisional 
release  (CH)  mechanism  of  dissociation,  defined 
by  .Ar-f-llro — ►.Ai'Br- b Hr,  contributes  to  overall 
di.s.soeiation  otdy  at  l.')lll)°K.  The  absence  of  the 
CK  mechanism  at  a higher  temiierature  is  cer- 
tainly due  to  the  weak  attraetion  between  argon 


and  bromine  atoms.  The  complete  collisional 
dissociation  (CCD)  mechanism,  defined  by 
Ar-+-Hr.2— >Ar-|-2Hr,  predominates  at  all  tempera- 
tures. 


5.  Energy  and  Angular  Momentum  Transfers 
in  Collisions 

The  average  energy,  transferred  per  collision, 
(A/i  ),  is  .shown  in  Fig.  .'i  as  a function  of  Hpi 
internal  energy.  In  a reactive  collision  with  .a 
bound  molecule,  (AE)  (reactive),  is  positive. 
Furthennore,  (AE>  (reactive)  is  significantly 
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Kio.  (').  Avonigc  luiRular  momi'iitiim  triiMsfer  ts.  reagent  internal  energy.  For  notations,  see  captions,  Fig.  ,5. 


larger  tlian  the  minimum  amount  of  energj- 
necessary  to  excite  the  moieeiile  just  u|)  to  the 
(li.ss<ieiation  limit.  However,  if  the  <liatomic 
molecule  i.s  in  .a  highly  metast.alile  state,  (AE) 
(reactive)  is  negative.  Since  de-exeitation  of  the 
vihr.ational  motion  alone  cannot  dissociate  a 
classically  stahle  molecule,  the  dissociation 
involves  loss  of  rotational  energy,  which  leads  to 
the  lowering  of  the  rotational  harrier. 

In  order  to  sluily  further  the  angular  mo- 
mentum tran.sfer  in  energize<l  Hr.,  molecules, 
we  defined  the  average  angular  momentum 
transferred  in  dis.soeiative  collisions  ns 

I t 

in  which  //  is  the  orhilal  angular  momentum 
of  tlie  di.ssociated  pair  after  collision,  i.s  the 


rotational  angtilar  momentum  of  the  initial 
diatom,  w;  is  the  appropriate  weight  of  the  fth 
trajectory,  and  5,''  is  a delta  function,  which  is 
unity  if  the  trajectory  is  dis.'ociative  and  zero 
if  otherwise.  The  average  angular  momentum 
transfer  in  non-dissociativc  collisions  is  antdo- 
gously  defined;  i.e. 

(Af  )nd=H  (//— 

I t 

where  1/  i.s  the  rotational  angular  momentum 
of  the  diatom  after  a non-dissoeiative  collision. 
Thus  defined,  (SI)  is  jiositive  if  there  is  an  in- 
crease of  aitguliir  momentum. 

The  values  of  (A/’  for  dissociative  as  well  as 
non-di.ssociative  collisions  were  plotteil  against 
the  His  internal  energy  tl'ig.  (i).  For  the  hound 
.states  (A/'d  was  po.sitivc,  indicating  rotational 
excitation  during  the  dissociation.  I'or  metastahle 
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Fir..  7.  (h)  FvolutioM  of  («')  the  average  internal  energy  in  the  reactive  energy  range,  tb)  Involution 
of  the  reactivity,  a-s  measured  by  (S ),  the  dia.soeiative  croas-seetion  iu  the  reactive  energy  range.  For  (>0(K)°K, 
the  ahsei.asa  is  contracted  by  a factor  of  two. 

species,  (Al)d  is  negative,  indicating  the  lower- 
ing of  tingular  momentum  during  dissocititive 
collisions.  This  finding  corroborates  the  (Ai’) 
plots  (Fig.  5),  and  again  indicates  that  rotational 
de-excitation  contributes  to  tlie  dissociation 
of  inetastaide  molecules.  Collisions  of  metastable 
molecules  with  internal  energy  below  -|-2kT, 
yield  negative  (M)d  but  positive  (SE>  beeau.se 
at  this  energy,  the  e.xeitation  of  vibrational 
motion  outweighs,  on  the  average,  the  de-ex- 
citation of  the  rotational  motion.  Figure  (>  sub- 
stantiates an  earlier  conclusion-^  that  vibra- 
tional-rotational coupling  is  strong  if  molecules, 
with  energy  close  to  the  dis.sociation  limit, 
collide  with  a third  body. 

For  non-di.s.soeiativc  collisions,  Fig.  0 indicates 
that  reagent  molecules  in  highly  excited  states 
also  experience  rotational  de-excitation.  Further- 
more, inela.stic  collisions  yield  but  a very  small 
(AE)„d  and  at  all  zt/RT  values,  while 

the  corresponding  (A£,’/d  and  {M)d  are  con- 
siderably larger  (Figs,  o and  6). 

6.  Non-Equilibrium  Effect  in  Dissociation 

The  non-equilibrium  effects,^* di.scu.s.sed 
previously,^**  occur  because  reacting  sys- 
tems, under  a steady  state  condition,  cannot  be 
uniquely  defined  by  a Boltzman  distribution  at 
temperature  T.  The  non-equilibrium  effects  may 
be  characterized**  by  the  ratio  of  the  rate  con- 


stant at  .steady  state,  i.e.  k/‘,  to  A'/'’  and  can  be 
expressed  as 

f=  k/‘  k,^= J ihi  j ^ 0">(«)  (hi 

■ (^) 
Two  factors  contribute  to  the  non-eiiuilibriuin 
effect.  In  the  first  place,  the  reactive  energy  range 
is  (lei)leted  by  reactive  collisions.  This  factor  is 
determined  by  the  ratio  of  the  integrals  of  the 
di.stribution  functions  [e<i.  (4)].  Furthermore, 
the  reactivity  in  the  reactive  zone  is  also  ile- 
crea.sed;  this  factor  is  determined  by  the  ratio  of 
average  cross-sect  ions. 

Our  method  of  calculation  of  / factors  is  based 
on  multiple  collision  aiiproach,**  which  proved  to 
be  successful  in  earlier  trajectory  studies  on  re- 
combination."’-*"*^ The  time  evolution  of  the 
])opulatioi'  in  the  reactive  zone  can  be  expressed 
as 

dAo,  0t=  — Jd—  /d-.nd+  Jnd-.d 

in  which  Jd  is  the  dissociation  flux,  and  Jd^^ 
and  are  fluxes  representing  de-excitation 

and  excitation  of  molecules  from  and  into  the 
reactive  energy  range.  In  the  above  expression, 
flux  due  to  recombination  of  atoms  is  neglected. 

In  the  present  calculation,  we  generated  an 
equilibrium  ensemble  of  diatomic  molecules  at 
temperature  T with  internal  energy  within  the 
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roMCtive  /.nil  ■.  \Vc  tlicli  siilijcctcil  tlii>  iiiitiiil 
t'lisi'iiilili'  tn  I'nllisinii.'  witli  inert  (;ns 

atniiis  111  tlie  sniiie  teiii|)ei;it me  T.  A>  ii  result  nf 
eollisiniis,  a fraetioii  of  tlie  ensemble  (li.ssneiateil, 
and  the  rest  underwent  inelastic  eiiei'tjy  transfer. 
We  a.s.suine<l  that  the  transfers  nf  ninleeiiles  in 
and  nut  nf  the  reaelive  reninii  resulting  fmni 
I inelastic  eulhsiniis  are  identical,  i.e. 

i ,/„j  „i.  This  assiiinptinii  implies  that  the  ilepletinn 

i of  the  reactive  /.one  is  solely  caused  by  dissoeia- 

‘ tion,  but  ehanttes  of  the  eneiyty  distribution 

within  the  reactive  /.one  are  the  re.siilt  of  inelastic 
collisions  as  well  as  dis.sociation.  .At  a steady 
state,  the  eiieifty  distribution  of  the  en.'enible, 
and  its  reactivity,  are  invariant  with  re.spect  to 
time  and  number  of  collisions.  Therefore,  the 
average  internal  energy,  in  the  reactive  /.one, 
(a'h  and  the  average  dis.sociation  cross-section, 
(S),  of  the  ensemble  were  chosen  as  the  two 
liarameters  in  monitoring  the  e.'tablishment  of 
the  steady  state. 

r.sing  this  model,  we  calculated  the  iion- 
eipiilibriuni  effects  in  dissociation  at  ihoOtl  and 
■ ()t)IIO°K.  .At  each  temperature  an  eight  hundred 

trajectory  sample  per  collision  was  used,  but 
stratified  and  importance  .sampling  was  not 
attempted.  The  average  internal  energy  and  the 
average  dissociation  cross-section  in  the  reactive 
energy  range  were  studied  as  a function  of  the 
number  of  collisions.  It  jiroved  difficult  to  deter- 
mine unambiguously  at  which  point  the  steady 
state  is  e.stabli.shed.  Kandom  fluctuations  were 
pronoimced,  particularly  in  (.S')  vs.  number  of 
collisions  plots.  Therefore,  the  ensemble  evolu- 
tion was  studied  uii  to  twenty-nine  collisions 
at  7’=3500°K  and  up  to  forty-six  collisions  at 
(iOOO°K.  It  was  possible  to  obtain  relatively 
smoothly  averaged  curves  (Fig.  7).  The  non- 
equilibrium correction  was  then  obtained  u^ing 
equation  (4). 


7.  Comparison  with  Experimental  Data 

The  calculated  equilibrium  dis.sociatioii  rate 
coefficients  yield  an  excellent  linear  .Arrhenius 
plot,  with  an  activation  energy,  A/','(act),  of  44. .A 
Kcal mole,  assuming  that  the  |)re-exponenlial 
factor  in  the  rate  constant  is  temperature  inde- 
pendent. This  number  is  within  one  Kcal  mole 
of  the  di.ssociation  energy  of  lii'2  (4.A..A  Kcal  mole). 
That  the  activation  energy  for  A',f'  is  approxi- 
mately eipial  to  the  dissociation  energy  of  the 
diatomic  molecule,  was  al.M)  reported  by  previous 
invest  igators.’'- 

If  the  non-etjiiilibriiim  corrections  are  made, 
the  computed  dissociation  rate  constants  show 
an  appreciable  negative  curvature  in  the  .Ar.- 


heniiis  plot.  Thi“  curvature  is  -oniewhat  more 
pronounced  at  higher  temperature-.  The  average 
slope  for  such  a plot  yield-  A/'.'(acl i=  hi  kcal 
mole. 

In  order  to  compare  the  pre-eni  calculations 
with  the  results  obtained  by  flash  photoly-i-"’ 
and  shoi'k  wave  techniques,' we  converted 
the  dissociation  rate  con-taiil.-  into  recombina- 
tion rate  I'onstaiits  using  appropriate  eipiilibrii'"! 
constants  (Fig.  3).  Since  the  non-eipiilibriiim 
correction  factors  were  calculated  only  at  fl.'ibb 
and  t)00l)°K,  it  was  necessary  to  estimate  Ihe-e 
factors  at  l.allO  and  l!.a(IO°K  by  extrapolation 
between  3,")00°K  and  lower  temperature  data.'-" 
Such  extrapolation  proved  to  be  iinambigioiis. 

It  follows  from  Fig.  2 that  the  equilibrium  re- 
combination rate  coefticients  yield  a positive 
teniiieratnre  dependence.  .After  non-equilibrium 
effects  were  included  in  calculations,  the  rei'om- 
bination  rate  constants  obtained  in  the  present 
calculations  were  found  to  be  consistent  with 
flash  photolysis  data.'" 

Since  the  calculations  of  non-eipiilibrinm 
correction  factors  were  done  only  at  two  tempera- 
tures, it  is  difficult  to  estimate  the  precision  of  the 
calculated  rate  constants.  It  is  iirobaliie.  however, 
that  the  rate  constants  are  reproducible,  ap- 
proximately within  a factor  of  two.  For  this 
rea.-on,  the  rate  con.staiits,  obtained  in  all  ‘‘reli- 
able" shock  wave  experiments,  are  also  consi.s- 
tent  with  calculated  dissociation  and  recom- 
bination data. 

Finally,  it  may  be  iiointed  out  that  the  agree- 
ment between  trajectory  calculations  in  recom- 
bination domain  (open  triangles)  with  the  tra- 
jectory calculations  in  the  dissociation  domain 
(open  .squares),  demonstrates  well  the  validity 
of  the  rate  quotient  law  [A(forward),''A'(reverse)  = 
A'(e(|l3  under  the  steady  state  condition  (Fig.  2). 
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COMMENTS 


IF.  /■'ors/,  r Hirer. 'till'  l.mni,  Canada.  1 wonder 
if  rate  eonstaiit.s  mea.sureil  in  various  shock  tube 
ami  Hash  [ihotolysis  e.yperimeuts  are  really  true 
skruly-sUile  rate  constants.  Could  it  be  that  the 
di.sciepancy  amonn  the  eN|)criinental  data  of 
varioirs  investigators,  alluded  to  in  your  paper,  is 
line  to  measurements  done  at  various  stages  of 
the  initial  transient  that  precedes  the  establish- 
ment of  the  steady-state? 

Aulhurtf  Rcpb).  In  order  to  measure  rate  cou- 
.stants  in  Hash  iihotolysis  experiments,  the  con- 
centration of  diatomic  molecules  is  followed  over 
a relatively  long  time.  In  these  experiments,  the 
rate  coellicient  is  always  a time-independent  rate 
constant.  Therefore,  there  seems  to  lie  no  reason 
to  suspect  that  in  these  ex|ieriments  the  .steady 
state  is  not  reaclied.  On  the  other  hand,  the  rate 
Ilf  dissociation  of  diatomic  molecules  is  olten 
determined  from  the  initial  .slo/ie  in  .sliock  tube 
experiments.  This  slope  may  be  determined  over 
ii  relatively  short  period  of  time  of  a few  inicro- 


seconds.  It  may  well  be  that  the  ili.se  re  pa  ncy 
between  some  shock  wave  data  and  Hash  photol- 
ysis data  is  due  to  the  inahilit.v  of  a .shock  wave 
experimentalist  to  measure  true  steady-state  rate 
constants.  Detailed  calculations'  demon.strate 
that  such  a supposition  would  indeed  reconcile, 
some  .shock  wave  and  Hash  photolysis  data. 


liEFEItKNCF, 

1,  le,  ,1.  K.  K.  AM)  Hvaxs,  G.:  ./.  Chem.  Phys.,  At. 
3425  (lOtill). 

4 

Jlajime  .\Iii/nmii,  Toraij  Industries,  Inc.,  Jaiian. 
Can  you  exiiect  much  more  contribution  of  colli- 
sioual  release  if  you  change  3rd  body  from  .\r  to 
other  species,  or  change  reactant  from  HPi  to  h. 
according  to  your  calculation? 


KLKMKNTAHY  H1;A(  TIONS 


Authors'  Reply.  Tho  collisiomil  release  mecha- 
nism is  more  important  if  dissociation-recombina- 
tion reactions  occur  at  lower  temperatures  and 
involve  stronger  interactions  between  the  recom- 
bining atom  and  a third  body.  If  the  third  body 
is  an  inert  gas  and  the  temperature  is  sulliciently 
high,  so  that  the  dissociation  reaction  is  ai>- 
preciable,  the  collisional  release  mechanism  is 
not  important.  For  the  case  of  iodine  dissociation 
this  point  was  demonstrated  recently.'  On  the 
other  hand,  at  room  temperature  the  atomic 
recombination  apijear;  to  proceed  via  a collisional 
release  mechanism,  provided  that  the  third  body 
is  sufficiently  heavy,  such  as  Xe.^ 


HKFFllKN’Ci;.s 

1.  WoNO,  W.  II.  AND  Hunss,  I'roc.  Itoy.  .Soc, 
bond.  (In  prcs.s). 

2.  WoN<i,  W.  II.  AND  Burns,  <■.:  J.  Oheni.  Phys. 
5H,  44.-.P  (197.1),  J.  Chem.  Phys.  6t>.  2974  (1971). 


C.  I‘.  Quinn,  Shell  Research  Ltd.,  England. 
Does  your  lasscr  schlieren  techtdipie  allow  you 
to  Isolate  an  incubation  |H-riod  for  bromine 
dissociation  and  if  .so  how  dues  if  compare  with 
the  vibrational  rela.xation  time? 

Authors’  Reply.  W'e  were  not  able  to  detect 
an  incubation  time  for  Br-j  in  our  exiK‘riments. 
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